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Abstract 

The harmonium model has long been regarded as an exactly solvable laboratory 
bench for quantum chemistry [1]. For studying correlation energy, only the ground 
state of the system has received consideration heretofore. This is a spin singlet state. 
In this work we exhaustively study the lowest excited (spin triplet) harmonium state, 
with the main purpose of revisiting the relation between entanglement measures and 
correlation energy for this quite different species. The task is made easier by working 
with Wigner quasiprobabilities on phase space. 

1 Introduction 

Replacing the wave function of electronic systems by the reduced 2-body density matrix 72 
tremendously saves computation without losing relevant physical information. Till very 
recently, the solutions to the iV-representability for that matrix [2,3] were impractical. This 
certainly did not impede great advances in the use of 72 for many-electron quantum systems 
— see for instance [4]. Now a constructive solution [5] to that representability problem, 
leading to a hierarchy of constraints [6] on the variation space for 72, has been unveiled. 



At any rate, the last fifteen years have witnessed a justifiable amount of work in trying 
to obtain the 2-body matrix as a functional of the 1-body density matrix 71. Starting with 
the pioneer work by Miiller [7] , several competing functionals have been designed, partly out 
of theoretical prejudice, partly with the aim of improving predictions for particular systems. 
We shall discuss pure state representability for 7x in the case of our interest in Section 6. 

Two-electron systems are special in that 72 is known "almost exactly" in terms of 7^ Let 
us express 71 by means of the spectral theorem in terms of its natural orbitals and occupation 
numbers. For instance, the ground state of the system admits a 1-density matrix: 

-yi(x,x') = (titi'+Wi'hi^O = (titr+W>i0^^&( r )#(O- (1) 

i 

Here J2i n i = 1- Mathematically this a mixed state. The corresponding 2-density matrix is 
given by 

l2 (x 1 ,x 2 ;x , 1 ,x , 2 ) = (til 2 -lit 2 )(ti4 2 ' -li't 2 O^^^(n)0 l (^ 2 )0*( r i)^( r 2) ) 

with coefficients q = ±y/nl. (2) 

The expression is exact, but the signs of the q need to be determined to find the ground 
state [8,9]. Note that 7! = 72. The first excited state of the system admits a reduced 

1- density matrix of the kind: 

li(x;x ) = (spin factor) x (<Mr)^(0 + &i+i( r )02i+i(O) 

ij 

with J2i n i = 1 an d spin £ {titi'j §(titi' + 4i4i')> liir }■ The corresponding spinless 

2- density matrix 72(7*1, r 2 ; r[, r' 2 ) is given by 

^^[0 2 *(ri)0 2m (r 2 )0^ 

ij 

- 02i(T*i)02i+i(r* 2 )02j+iW)02j(^2) - 02i+i(n)0 2 i(r 2 )0; i (r / 1 )0; j+1 (r;)], 
with coefficients q = +y/ni. 

Due to the antisymmetry of this state, there is no ambiguity in the choice of sign. 

A completely integrable analogue of a two-electron atom, here called harmonium, de- 
scribes two fermions interacting with an external harmonic potential and repelling each 
other by a Hooke-type force; thus the harmonium Hamiltonian in Hartree-like units is 

tt Pi 1 V\ . k, 2 , 2\ ^ 2 /o\ 

H = y + y + +r 2) ~ 4 r i2- (3) 

where r\ 2 := \v\ — r 2 \. This model is rooted in the history of quantum mechanics: Heisenberg 
first invoked it to approach the spectrum of helium [1] . 
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Several problems related with this model — although not quite the present one — are 
analytically solved; and so it is tempting to employ it as a testing ground for methods used 
in other systems, such as the helium series. Indeed, Moshinsky [10] reintroduced it with 
the purpose of calibrating correlation energy. There is considerable interest nowadays on 
learning from harmonium, including further study of correlation [11-13], approximation of 
functionals [14,15], and beyond quantum chemistry, questions of entanglement [16-19] and 
black hole entropy [20]. 

In the past, harmonium problems have been attacked with ordinary wave mechanics [21]. 
Now, for the analysis of harmonium the phase space representation of quantum mechanics 
recommends itself. The deep reason for this is the metaplectic invariance of that formalism 
[22], hidden in the standard approach: this made it easy to solve the sign dilemma in the exact 
Lowdin-Shull-Kutzelnigg formula [8,9] for 72 in terms of 71, for two-electron systems [23,24]. 
We come to this at the end of the next section. Such a phase-space description was taken up 
first by Dahl [25], and then developed, within the context of a phase-space density functional 
theory (WDFT), by Blanchard, Ebrahimi-Fard and ourselves [23,24,26-28]. 

Our goal in this article is to understand, in WDFT terms, the first excited state of 
harmonium. As for helium-like atoms, we expect it to be the lowest spin triplet state, 
to which we refer simply as the triplet. Particularly we make clear the nonexistence of a 
phase dilemma in this situation, and pinpoint the similarities and differences between the 
relative behavior of entropy and correlation energy for the (spin singlet) ground state and for 
the triplet. Again, and essentially for the same reason, WDFT shows its worth here — see 
Section 6. 

The customary plan of the paper follows. In Section 2 we briefly recall for the benefit of 
the reader our treatment for the singlet ground state; this helps to introduce the notation. 
Sections 3 and 4 deal with the general mathematical structure of triplet 1-body Wigner 
functions. Section 5 computes the Wigner quasiprobabilities for the harmonium triplet. 
Section 6 deals with the corresponding natural orbitals. In Section 7 the behaviour of 
the occupation numbers, obtained numerically, is compared to that of the ground state. 
Section 8 continues this comparison in the setting of quantum information theory. The 
relative correlation energy for the triplet is smaller than for the singlet, just as is the purity 
parameter. The proportionality between entropy and correlation energy, observed in the 
weak correlation limit for the singlet, fails for the triplet state. Section 9 is the conclusion. 

2 Wigner natural orbitals for the harmonium ground state 

Given any interference operator |^)($| acting on the Hilbert space of a two-electron system, 
we denote 



p2*$(ri, r 2 ; Pi,P2, ft, ft; ft', ft') 

:= I *(n - zi, r 2 - z 2 ; q u q 2 ) + z u r 2 + z 2 ; e a(pi-*i+w«) dz 1 dz 2 . 




(4) 
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These are 4x4 matrices on spin space. When \P = $ we speak of Wigner quasiprobabilities, 
which are always real, and we write d 2 for P 2 . The extension of this definition to mixed 
states is immediate. The corresponding reduced 1-body functions are found by 

P\^{r\\V\S\\<;v) = 2 J P 2 **(ri, r 2 ; Pi,p 2 ; ft) dr 2 dp 2 dq 2 . 

These are 2x2 matrices on spin space. When \1/ = $ we write c?i for Pi. The associated 
spinless quantities d 2 (ri, ^2; P11P2) and di(r;p) are obtained by tracing on the spin variables. 
The marginals of d 2 give the pairs densities p 2 (ri, r 2 ), 7r 2 (pi,p 2 ). The marginals of di 
give the electronic density, namely p(r*i) = f di(T*i,Pi) dpi, and the momentum density 
7r(pi) = / di(ri,pi) dri. 

It should be obvious how to extend the definitions to iV-electron systems and their reduced 
quantities; the combinatorial factor for djy >->■ d n is (^). 

Putting together (2) and (1) with (4), one arrives [23] at: 

d2(ri,r 2 ;pi,p 2 ;<ri,<r 2 ;<^,<r 2 /) = (spin factor) x ^ ^ Xij(n;Pi)Xij(T- 2 ;p 2 ), (5) 

and di(ri;pi;^i,^/) = 2 ^ d 2 (ri, r 2 ;p 1 ,p 2 ; q 1 , q 2 ; q v , q 2 ) dq 2 dr 2 dp 2 

= (titr + iiiv) Ui Xi(ri;pi). 

i 

Here rii are the occupation numbers with ^ f rij = 1, the Xij the natural Wigner interferences 
and Xi := Xm denote the natural Wigner orbitals; the spin factor is that of (2). Evidently 
(titr + 4i4i') is a rotational scalar. We replace it by 2 in what follows. 

The relation q = holds. In principle there still remains the problem of determining 

the signs of the infinite set of square roots, to find the ground state. To recover d 2 from d\ is 
no mean feat, since it involves going from a statistical mixture to a pure state — see below. 

Bringing in extracule and intracule coordinates, respectively given by 

1 1 

R=-j={r 1 + r 2 h r = —(r 1 -r 2 ), 

P = ^=(Pi +P2), P=^(Pi-P 2 ), 

the harmonium Hamiltonian is rewritten: 

u u u P 2 u 2 R 2 p 2 fi 2 r 2 

H = H R + H r := 1 h — + - — . 

R-r r 2 2 2 2 

We have introduced the frequencies u := y/k and /i := \Jk — 5. Assume 5 < k, so both 
"electrons" remain in the potential well. For the harmonium ground state the (spinless) 
Wigner 2-body quasiprobability is readily found [25]: 

d 2 {r u r 2 ,p u p 2 ) = ^expl —I expl — —I. (6) 
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The reduced 1-body phase space quasiprobability for the ground state is thus obtained: 
dl (r 1 ;p 1 ) = -( ^^y /2 e - 2 ^/^e- 2 ^/^). 

■K 6 \(C0 + fl) 2 J 

For its natural orbital expansion, with i integer > and Li the corresponding Laguerre 
polynomial, one finds [23] 



4 = n,= , ( ) 21 =■■ (1 - t 2 ) f » ; (7) 



( v ^ +v ^) 2 Vv^+v^ 



= fi(xi,Pix)fi(yi;pi y )fi(zi;piz), where 

fi(x; Px ) = -(-irL^2^a: 2 + 2py^7i) e ^- 2 ^/^. 

7T 

The functions determine up to a phase the interferences: for j > k, 



1 - -NfcvA;! , 2 2 .(j- k )/2 



f jk (x, Px ) = -(-l) fe ^= (2 v ^7i a ; 2 + 2p 2 / v ^7i) 1 

x e~ i{j - k)l£> Li~ k (2y/UJjlx 2 + 2pl/y/LJj2)e-^ x2 - p " / ^, 

where 

i? := arctan^/^/u^Ltx). 

The are associated Laguerre polynomials. The are complex conjugates of the f jk . 
Now, with the alternating choice (unique up to a global sign): 

Cl = (-) l ^ = V\^(-t)\ 



and the above fjk, formula (5) does reproduce (6). This was originally proved in [23], and 
verified by minimization in [24] ; we refer the reader to those papers. Trivially, the same sign 
rule holds for natural orbitals of the garden variety (2). 



3 Generalities on the triplet state 

For a general two-electron system in a triplet spin state the reduced 1-density possesses three 
different spin factors, say 

titr and |(titi' +4-1^1') and ^ v . 

While the spatial function for the ground state is symmetric, and consequently its spin part 
antisymmetric, for the first excited state the situation is exactly the opposite: the spatial 
function is antisymmetric and its spin part is symmetric. This leads to important differences 
between both cases for the natural orbital decomposition. 
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General triplet states are describable in the form [8,21]: 

*ti(ri,r 2 ;ft,fe) = tit 2 Yl \ Ci i [^(n)^(^) - ^j(ri)^i(r 2 )], 

ij 

*to(»*i,r2;ft,fe) = -7= (T1I2 + lit 2 ) 2 Cij [^( r i)^'( r 2) -^■(''i^M, 
*t,-i( r i> r 2;ft>&) = lil 2 $^ ^Cij [^(ri)^(r 2 ) -^(ri)^(r 2 )], 



where Qj = — Cjj. Here {^i} is a complete orthonormal set. In the absence of magnetic fields, 
the wave functions can be taken real. We thus assume that the matrix C = [c^] is real, as 
well as the functions ipi. Wave function normalization gives rise to Tr(C* C) = c % — 1- 

For the spin part, a less conventional and more cogent description is found in terms of 
polarization vectors and the correlation tensor [29, App. F]; however, it is hardly worthwhile 
to introduce it here. So we shall be content with presenting the Wigner 2-body quasiproba- 
bilities for triplet states under the matrix form 



A^ti^ti — tit 2 ti't 2 ' d>2 



(d 2 0\ 





\0 0/ 



j ^2* t _i*t_i — iii 2 4i'i 2 ' d 2 



/o 0\ 





\0 d 2 J 



-P2*io* tt) — 2 (til 2 + iit 2 ) (ti42' + li't 2 ') ^2 = 2 



fo 0\ 

d 2 d 2 

d 2 d 2 

\0 0/ 



where d 2 is the spinless Wigner 2-body quasiprobability, given by the expression 
d 2 (r 1 ,r 2 ;p 1 ,p 2 ) 

= \ Yl Ci i Ckl J ~ Zl )^'( r2 ~ Z2 ) - i>i{r\ - zi)^i(r 2 - z 2 )\ 

ij,kl 

x [^*(n + z 1 )^(r 2 + z 2 ) - ^(n + zi)^(r 2 + z 2 )] e 2 ^ 1 '* 1 ^ 2 ' 22 ) d Zl dz 2 
= ^^CijCki [P ik (r 1 ]p 1 )P j i(r 2 ]p 2 ) - P u (r 1 ;p 1 )P jk (r 2 ;p 2 ) 

ij,kl 

- P jk (ri;pi)Pii(r 2 ;p 2 ) + P ji (ri;pi)P ife (r 2 ;p2)]- 
By integrating out one set of coordinates, we obtain the 1-body quasiprobabilities: 

Pl^n^n = tf di = ^q 1 , -Plt t ,-i*t,-i = II' di = j^J , 

Pi * t0 * t0 = \ (tf + II') di = ^ ( Q 1 • 



(8) 
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Here d\ is the spinless 1-body quasidensity corresponding to the triplet: 
di(r;p) = 2 J d 2 {r,r 2 ;p,P2) dr 2 dp 2 

= \ ^CijCki J [P ik {r;p)P j i(r 2 ;p 2 ) - Pu(r;p)P jk (r 2 ;p 2 ) 

ij,kl 

~ Pjk(r;p)Pu(r 2 ;p 2 ) + Pji(r;p)P ik (r 2 ;p 2 )] dr 2 dp 2 
= \Y1 Ci i c ki[ p ik(r; p) 6{ - Pa(r; p) 5{ - P jk (r; p) 5\ + P^r; p) 5\] 

ij,kl 

= 2^2a k c jk Pij(r; p) =2^2 d ij p ij( r ^ P)» 

ij,k ij 

where D = CC l = —C 2 is a positive definite matrix. 

4 The Schmidt decomposition of the triplet 

Let C be any real antisymmetric square matrix. It is well known that there exists a real 
orthogonal matrix Q such that A = Q l CQ, with A a real block-diagonal matrix: 

A = diag[A ,A 1 ,...], A) = 0, A *=(-a i o)' 
By convention, here Oj > 0. Therefore 
^CijCMP^r^p^Pjiir^p^ 

ij,kl 

[Qi,2vQj,2v+l — Qi,2v+lQj,2v\ [Qk,2wQl,2w+l ~ Qk,2w+lQl,2w] Pik{ri;Pi)Pji{r 2 ;p 2 ) 

ij,kl,vw 

= ^2 a va w [qi, 2 vPik(ri;pi)q k , 2w qj, 2 v+iPji(r 2 ;p 2 )qi t2w+ i 

ij,kl,vw n f \ id ( \ 

- qi,2vPik{ri;Pi)q k , 2w+ iqj, 2v +iPji{r 2 ;p 2 )q lt2w 

- qi,2v+iPik(ri;pi)q k , 2w qj, 2v Pji(r 2 ;p 2 )qi, 2w+ i 
+ qi,2v+iPik{ri;Pi)q k , 2w+ iqj, 2 vPji{r 2 ;p 2 )qi,2w\- 

Let us now make the definition Xr P (r;p) := J2mk<irnr P m k(r;p) q kp , so that P mk (r;p) = 
^rpQmr Xrp(r; p) qkp- This is the set of Wigner natural orbitals, and has the following nice 
property: 

J Xr P (r; p)drdp = J ^ qm r P m k(r; p)q kp drdp = ^ Qm r qk P 8™ = 5 r p . 

mk mk 
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Hence, 

^2c ij c kl P ik (r 1 ;p 1 )P j i(r 2 ;p2) 

ij,kl 

= ^2 a v a w [X2v,2w(r 1 ;p 1 )x2v+i,2 w +i(r 2 ; p 2 ) - X2v,2w+i(ri,Pi)x2v+i,2w{r 2 ;P2) 

WW 

- X2v+i,2w(r 1 ;pi)x2v,2w+i(r 2 ; P2) + X2v+i,2w+i{r 1 ; Pi}X2v,2 W {r 2 ] p 2 )] . 
The other three summands in (8) yield the same expression. For instance, the third is 

- c ij c ki p u( r ^ Pi) p jk(r 2 ; p 2 ) 

ij,kl 

= — ^ 0, v a w [qi,2vQj,2v+l — Qi,2v+lQj,2v] [Qk,2wQl,2w+l — Qk,2w+lQl,2w] Pu(fi', Pi) Pjk(r 2 ; p 2 ) 

ij,kl,vw 

= - ^ avaw[qi,2vPu(ri;pi)qi,2w+iqj,2v+iPjk{r 2 ;p 2 )qk,2 W 

ij,kl,vw 

- qi,2vPii{ri;Pi)qi, 2w qj,2v+iPjk{r 2 ] p 2 )qk,2w+i 

- qi,2v+iPu(ri;pi)qi, 2w +iqj, 2v Pjk(r 2 ;p 2 )qk,2w 
+ qi,2v+iPii{ri;Pi)qi, 2w qj, 2v Pjk{r 2 ;p 2 )qk,2w+i] ■ 

This leads to the same contribution as the first summand. Then use symmetry under the 
interchange of the two particles. In summary, 

d 2 {r u r 2 ;p^p 2 ) 

= ^ a v a w [X2v,2w{ri;pi)x2v+i,2w+i{r2; P2) - X2v,2w+i{ri,Pi)X2v+i,2w{r 2 ;p 2 ) 

ww 

- X2v+i,2w(ri;Pi)x2v,2w+i(r2]P2) + X2v+i,2w+i(r 1 ;p 1 )x 2v ,2w(r 2 ;p 2 )]. (9) 

The reduced 1-body phase space (spinless) quasidensity for the triplet is obtained, as 
before, 

d 1 (r;p)=2 / d 2 (r, r 2 ;p,p 2 ) dr 2 dp 2 = 2 ^ a 2 w [x(2w,2w)(r;p) + X(2w+i,2w+i)(r;p)\. (10) 

J w 

Notice that in the previous equation each occupation number rii := 2af appears twice. This 
is a consequence of the Pauli exclusion principle. 

Unlike the singlet case, there is no sign rule to be deciphered here. Instead there are the 
ambiguities: 

X2w,2w = X2w,2w COs2 @w — (X2w,2w+1 + X2w+l,2w) S ^ n @w COS 9 W + X2w+l,2w+l Sm ®wi 
X2w+l,2w+l — X2w,2w S i n2 ®w + {X 2w ,2w+l + X2w+l,2w) COS 9 W + X 2w +l,2w+l COs2 @w 

They clearly leave the form (10) untouched. We see here the action of SO (2) on each 
invariant block. One may choose the angles as to maximize their overlap with the leading 



natural orbitals for the ground state, as done in the seminal paper by Lowdin and Shull [8]. 
We omit that. Let us define 

A (cos9 w -sinftA 
w ' I sin 9 W cos 9 W j ' 



The above transformation can be construed as 



/ cos 9 V cos 9 W — cos 9 V sin ft 



X = (Ay ^ Ay,) X ' = 



i v cos 
cos 6 V sin 9. 



w 



cos 9 V cos 9, 



sin 9 V cos 9 W — sin 9 V sin 9, 
\ sin 9 V sin 9 W sin ft, cos ft 



w 



— sin 9 V cos ft 

— sin 9 V sin ft 
cos 9 V cos ftu, 
cos sin Q w 



sin ft, sin 9 W \ 

— sin 6^ cos 9 W 

— cos ^t, sin 9 W 
cos cos 9 W J 



x', 



with 



x ■-- 



( X2v,2w \ 
X2v,2w+l 
X2v+l,2w 
VX2u+l,2tu+iy 



and similarly for x', 



in the case t> = w. 

To similarly examine the symmetry of expression (9), again one does not have to contend 
with the whole tensor product matrix, since most contributions vanish. As regards the sum 
in (9), one can write in compressed form: 



/ cos 2 9 V cos 2 ft 

2 n • 2 . 



XX 



— cos 2 9 V sin 2 ft 
~ : - 2 " cos 2 ft 



- sin 2 ft 



I 



\ sin 2 9 V sin 2 6> w 



— cos 2 ft, sin 2 # w 
cos 2 9 V cos 2 ft^, 
sin 2 fty sin 2 ft,, 

— sin 2 9 V cos 2 6* w 



— sin 2 9 V cos 2 
sin 2 ft, sin 2 9 W 
cos 2 ft, cos 2 9 W 

— cos 2 0„ sin 2 9 W 



sin 2 fty 



,sin 2 ft„ \ 

— sin 2 ft, cos 2 ftt, 

— cos 2 9 V sin 2 
cos 2 9 V cos 2 ftj, / 



X X ) 



with 



/ X2v,2w(ri,Pl) X2v+l,2w+l(r2] P2)\ 

X2v,2w+i(ri;pi) X2v+i,2w(r2] P2 

X2v+l,2w(ri,p 1 ) X2v,2w+i(r 2 ; P2 

\x2v+i,2w+i(ri,Pi) X2v,2w(r2] P2) y 
One verifies that (9) is invariant under this set of transformations 



and similarly for XX- 



5 Lowest triplet state of harmonium 

The energy spectrum for harmonium is obviously (N + |)w + (N + |)//. Since ji < u, the 
energy of the first excited states is Ef s = (3u + 5/x)/2. For our present purposes, it is 
enough to choose an intracule excitation state along the a;- axis (say). The corresponding 
2-quasidensity is given by: 



2 / 2H R \ ( 2H r \f(p 1I -p 2x ) 2 + ^ 2 (x 2 -x 2 ) 2 1 
? exp(- — )<**(-— )( " ~ 2 
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Henceforth we work in the chosen nontrivial mode, since the problem factorizes com- 
pletely. By integrating one set of variables, the reduced one-body spinless quasidensity is 
obtained, after some work: 

3 



di(r;p) 



d 2 (r, r 2 ;p,p 2 ) dr 2 dp 2 



2 (2JuJJi 



r 2 2_ 2 

g Id + fl UJ + fJ,^ 



ur 2 + 



1 



71 + fi 

The marginals of d\ give the electronic density and momentum density: 

3 



-P' 



00 



(12) 



p(r) = / d 1 (r;p)dp 



2 2Joo~fi 



7T y oo + fi 

2 (2JoJfi 



^(p) = / d±(r; p) dr 



7T \00 + fi 

2 (2Jo~Jfi 



e "+m 



2^M r 2 



^ 2 



1 



7r(cu + /i) 



7T \0O + /i 

2 (2Joofl 



7T \0O + fi 



g " + M 
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2 

2"M , 

g 



1/2 



p 2 I dp 



oo 



2 . W + ^ 

cur + 



cur 2 + 



cu 



4cu 



rfr 



2 / 7r(cU + fi) \ ^ / CU + fi 1 



2cu/i 



Afi 



cu 



Finally, as expected, we get 

f 2 

7r(p) dp = p(r) dr 



71 \00 + fi 



7r(cU + /i) ^ X ^ 2 ^ 7r(cU + fi)^ 1 ^ 2 f CO + fi ^ 00 + fi 



2ufi 



Afi 



Aoo 



2. 



From the viewpoint of WDFT, the most interesting part of the energy corresponds to 
the interelectronic repulsion of this first excited state E 2 f s . The 1-body Hamiltonian is given 
by h(r,p) = p 2 /2 + cu 2 r 2 /2. It is a simple exercise to obtain the 1-body energy Ek s by 
integrating expression (12) with this observable: 

Elis ~2 + Afi ■ 

The interelectronic potential in (3) is (fi 2 — cu 2 )r 2 2 /4, so to obtain the repulsion energy E 2 { s , 
one has just to integrate expression (11) with this observable: 



R 



2fs 



[ 2 / 2H R \ 
J ^ eXP (-— J 



71 



2 exp 

2 , ,2 



cu ) / exp 





2H r \ 


exp^ 










fi ) 


2H r 


) 











2H r 
fi 



+ fir 



1 fi 2 — cu 2 
r\ 2 dR dr dP dp 



r 
~2 



drdp=\ — , 

4 fi 



which is 3 times the interelectronic repulsion energy for the corresponding mode of the 
singlet [23]. This is not surprising, since, in the triplet configuration the electrons tend to 
be mutually farther apart than in the singlet. 1 



1 Interestingly, (12) is a non-Gaussian Wigner function taking only positive values. This prompts two 
remarks. First, in consonance with common wisdom [30,31], it is confirmed that as of itself d\ is a nearly 
classical state. Second, there are mathematical recipes that produce such positive-valued Wigner functions 
representing mixed states [32]. It would be good to know whether or not (12) can be obtained as such an 
output. 
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6 Spectral analysis of the 1-body triplet state 

In order to determine the occupation numbers of this system, first we have to find the good 
coordinates. Let us perform the transformation 

(Q,P) '■= ((u;/i) 1/ ' 4 r, (w/i)" 1 ^ 4 ^) ; or, in shorthand, U = Su, 

where S is symplectic and u = (r,p). We may also write d := arctan(P/Q), so that 

P = Usind and Q = U costf. 

Recalling 2^/uJjI/(u) + /i) = (1 — t 2 ) / (1 +t 2 ) from (7), the 1-quasidensity (12) takes the simple 
form: 



7r(l + t 2 ) 3 
2(1 - t 2 ) 3 



e -(l-^)[/V(l+t 2 ) [/ 2 



1 + t 2 o 1 ~t ■ 2 o 

cos v H sin v 

1-t 1 + t 

1 + t 2 2t 



n(l + t 2 ) 3 

The one-body quasidensity may be expanded as follows: 



+ 



1-t 2 1-t 2 



cos 2^ 



'eft?. 



d^U, ti) = frs(U, ti) d rs where d rs = 2n I d^U, $)f* s (U, -&)U dU . 

rs 

Then, for r > s, 

2tt J f* s (U, ■&) di(C/, t?) C/ dZ7 dtf 

= *(1 + f 2 )3 ( ~ ir ^ g" (1 " t2)t/2/(1+t2) g" t/2 (2^ 2 ) (r ~ g)/2 ^r (2£/ 2 ) £/ 3 dt/ 



7T 



e 

— 7T 

2\3 



i(r— s)i? 



1 + t 2 



2t 



cos 2-^ 



1-t 2 1-t 2 

= 4 ( 1 ~ t2 ) J r e -(l-t 2 )Uy{l+t*) e -U* {2U 2 ){ r-s )/2 Lr(2C /2 )[/ 3 d[/ 

vr(l + t 2 ) 3 7 
"2(1 +t 2 ) 



X 7T 



i-t 2 r i-t 



so that 



d^U^) = J2dss(t)f ss (U,^ + d s+2 , s (t)f s+2>s (U^) + d s , s+2 (t)U s+2 (U^), 



where actually d s+2 ,s = d StS+2 . 

Using the standard Mellin transform [33,34]: 



li 



we obtain by fairly easy manipulations, 



d ss (t) = (l-t 2 ) 2 (st 2s - 2 + (l + s)t 2s ); 
d s , s+2 (t) = (1 - t 2 )V(s + l)(s + 2)t 2s+1 . 

This means that, to find the occupation numbers, one has to diagonalize a symmetric pen- 
tadiagonal matrix: 



2\2 



D = (1-t 2 ) 



( 1 





















1 + 2t 2 





ait 3 















2t 2 + 3t 4 





a 2 t 5 







; o 


ait 3 





3t 4 + At 6 





a 3 t 7 




' 





a 2 t 5 





At 6 + 5t 8 
















a 3 t 7 





5t 8 + 6t w 




V ; 















(13) 



where a s := a/(s + l)(s + 2) . 

It is readily checked that the trace of this matrix is 2, as it should be. Its eigenspaces split 
into two parts: £ 2 — Vi © V 2 , where V\ — { x : all x 2n = } and V 2 = { x : all a^n+i = }. 
They correspond respectively to the matrices 



D 



2\2 



= (1 - t 2 ) 



and 



odd 



2\2 



V 



/ 1 


aot 











A 


aot 


2t 2 + 3t 4 


a 2 t 5 













a 2 t 5 


At 6 + 5t s 


a 4 t 9 







o 





a^t 9 


6t 10 + 7t 12 


a 6 t 13 




i 








a 6 t 13 


8t 14 + 9t 16 • 




V ; 










7 


+ 2t 2 


ait 3 













c^t 3 


3t 4 + 4t 6 


a 3 t 7 













a 3 t 7 


5t 8 + 6t 10 


a^t 11 













a 5 t u 


7t 12 + 8t 14 


a 7 t 15 













a 7 t 15 


9t 16 + 10t 18 





\ 



It is easily checked that these matrices have the same set of eigenvalues, as they should, since 
the occupation numbers must appear twice. 

As was shown in Section 3, there is a skewsymmetric matrix C such that D = C l C. This 
matrix is tridiagonal, and is the sum of two skew-symmetric matrices whose diagonalization 
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is trivial: 



C 



/o 


-l 








• 




1 











• 













-V3t 2 


• 















• 
















• 




V 










■■) 



+ (l-t 2 



/o 











• 










V2t 





• 







-y/2t 








• 
















Vit 3 ■ 















• 




V 










■■/ 



=: A + B. 

Also, D is the sum of two Hermitian matrices, namely A 1 A + B t B, which is diagonal, and 
A l B + B l A. 

One is reminded here of the Weyl problem: given two n x n Hermitian matrices A, B 
whose spectra are known, what could the spectrum of their sum C := A + B be? Some facts 
are clear: with an obvious notation for the eigenvalues, these must satisfy 

ci H h c n = ai H \- a n + b 1 + h b n ; c\ < a± + b\\ 



less clear, but also true, are 



C2 < 0,1 + b 2 ; c 2 < a 2 + b\\ 



and so on. The conditions written above are already optimal for n = 2. The necessary 
constraints are all linear homogeneous inequalities, bounding convex polyhedra. Horn made 
a conjecture for the general form of such inequalities, which was eventually proved [35]. 

The pure-state iV-representability problem in quantum chemistry (or "quantum marginal 
problem", in the jargon of information theory) should be considered as solved, after the work 
by Klyachko [36,37]. It is of the same type and answered by similar inequalities. Both ques- 
tions reduce to finding moment polyhedra for coadjoint orbits of unitary groups (associated 
to pertinent Hilbert spaces), which are computed by Duistermaat-Heckman measures [38]. 
A very readable and up-to-date account of all this is [39]. The Hilbert spaces considered 
are finite-dimensional. However, the results are valid for finite-rank approximations in the 
chemical context, and the patterns of the inequalities extend in a rather obvious way. Thus 
it is scarcely surprising that the Weyl problem surfaces in this simple instance. We leave for 
the future consideration of the moment polytopes for the occupation numbers, 2 and choose 
in this paper a direct approach to the eigenpair problem, completed by numerical analysis. 

The matrices -Devon and D odd are tridiagonal symmetric real matrices. The general eigen- 
value problem for a matrix T of this kind reduces to solving the following set of recurrence 



2 The number of their extremal edges grows very quickly with N and the rank; this makes for precision, 
but also for strenuous work. 
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equations: 



(d Q h 

t\ d\ t 2 

t 2 d 2 h 

t 3 d 3 

u 

\\ ; ; ; 



/ O (n r )\ 

01 (n r ) 

02 (n r ) 

03 (n r ) 

04 (n r ) 



V 



/ 0o(n r )\ 

01 (n r ) 

02 (n r ) 

03 K) 

04 (n r ) 



/ 



V 



/ 



where n r is an eigenvalue. The general solution is completely given in terms of the occupation 
numbers, by the following formula [40, Sect. 5.48]: 



0m(A) 



0o(A) 

tit 2 ■ ■ ■ t r 



det[AJ — T] mm , for each m > 1, 



where [XI — T] mm is the upper left m x m submatrix of (A/ — T), and 0o(A) ^ is chosen 
so as to normalize the eigenvector. 

This result implies that T = QDQ\ where dij = n,^ is the diagonal matrix whose entries 
are the eigenvalues and q,ij = 4>i(rij). Since QQ f = Q l Q = 1, the following orthogonality 
relations hold: 



0m(n r )0;(n r ) = <J, m , m (n r )0 m (n s ) = 

r=0 m=0 

oo 

n r m (n r )0i(n r ) = rf m 5[" + £ m <^ m_1 . 

r=0 

In summary, for d\ we obtain 

- oo oo 
dl(-) = ^2 n r f2i,2i(-) 0cvcn,iK) + ^(/2i,2i+2 + /2i+2,2i)(-) 0cvcn,i (n r )0cvcn,i+l ( 



•i=0 



i=0 

oo 



+ ^2 hi+l,2i+l(-) 0odd,i( n r) + y^(/2i+l,2i+3 + /2i+3,2i+l) (') 0odd,j(n r )0odd,i+l ( n r) 



i=0 



Here n r depends solely on the parameter t of (7). 



7 Numerical analysis of the occupation numbers 

As advertised, to find the n r we fall back on numerical computation. Figure 1 shows the 
behavior of the rank-eight approximation of the eigenvalues, as t is varied. Note that the 
first eigenvalue is very close to 1 in the neighborhood of t — 0, while the others are very 
small. As the value of t rises, the first eigenvalue begins to decrease and the others rise for 
a while. In the neighborhood of t — 1 all eigenvalues approach zero. 

Note that t is a very nonlinear parameter: although t ~ 5/8 k for small 5, the value 
t — 1/2 means ji/u = 1/9 or 5/k = 80/81. This shows that, unless 5 is pretty close to the 
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Figure 1: First six eigenvalues of the matrix ^even- 
dissociation value, the harmonium triplet is not badly described by a Hartree-Fock state. 
Whenever t < 0.6, that is, S/k < 255/256, the first two occupation numbers contain almost 
all the physical information for the system. 

Also, one we can show that whenever t < 0.5, a good approximation to the five first 
occupation numbers is 

Ai « 1 - 3t 4 + 8t 6 , A 2 « 3t 4 - 8t 6 , A 3 « 5t 8 , A 4 « 7t 12 and A 5 « 9t 16 . 




Figure 2: First and second occupation numbers of the ground state and of the first excited 
state. 

Figure 2 compares the behavior of the first two eigenvalues for the singlet and triplet states 
of harmonium. In this sense, the Hartree-Fock approximation works better in the triplet 
case than for the singlet. Around t = 0.4 the second approximated occupation number for 
the latter is above 0.13, and for the former is below 0.052. The same behaviour was also 
observed in the toy model studied in [42]. This does not mean, however, that correlation is 
always weaker in the triplet state — see the next section. 
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8 Spatial entropy and correlation energies 



We move towards the comparison of the triplet system with the singlet system in regard to 
disorder (suppressing the spin variables). To measure this, a useful quantity is the linear 
entropy s associated to the 1-body function: 

s = i - n(di), 

where n(rf 1 ) is the purity of the system — see below. 3 Mathematically the quantity s is 
a lower bound for the Jaynes entropy, which has been used to quantify the entanglement 
between one particle and the other N — 1 particles of the system [42], and proposed as a 
handle on the correlation energies [43]. In this paper the singlet has been modelled in such 
a way that, for each one-dimensional mode: 

n gs ,i(di) = / d 2 1 (r;p)drdp = ^2n 2 l . 
J i 

Instead, for the triplet one should take for the excited mode: 

iW*) = i/*(r.;»>*,«h = E«?- 

This second definition is natural in that correlations due solely to the antisymmetric character 
of the wave function do not contribute to the entanglement of the system [18,44,45]. This 
ensures that the entropy for a 1-body function of the Hartree-Fock type is zero. 

In the singlet case, the occupation numbers are equal to (1 — t 2 )t 2% . Thus, the purity 
of this system is easily computable, to wit, n gs l (d 1 ) = (1 — t 2 )/{l + t 2 ) for each mode. 
This quantity coincides with the quotient of the geometric and arithmetic means of the 
frequencies, that is, n gs l = 2yJUJji/{uj + /i). For n modes one just takes the nth power [17]. 
Moreover, for small values of the coupling 8, we obtain 

- -~ --, (14) 



^ 32 w 4 ' 

which for this approximation is exactly the absolute value of the (dimensionless) correlation 
energy [26]. This appears to vindicate the contention of [43]. (Actually, for the singlet it is 
not difficult to compute the Jaynes entropy, given by 

Ei /-, 2x t 2 logt 2 
UilogUi = -log(l -t ) - . 

i 

This was done by Srednicki [20] some time ago.) 



3 Truth to be told, the notion of entropy native to the Wigner quasiprobability approach is the one 
discussed in [41]. We put aside the question of its eventual usefulness here. 
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For the triplet state, we have to compute Tr(df) for the matrix given in (13). Since 

/ l + a%t 2 a t(l + 2t 2 + 3t 4 ) • • \ 

(1 + 2t 2 ) 2 + a 2 t 6 

a t{l + 2t 2 + 3t 4 ) a 2 t 2 + {2t 2 + 3t 4 ) 2 + a 2 t w •• 



d\ = (1 - t 2 ) 4 
we get 



V 



7 



Tr(d 2 ) = (1 - t 



2\4 



= 2(l-t 2 ) 4 E 



4^a 2 t 2 ( 2i+1 ) + 2^ ? 2 t 4 ^ 1 ) 

8=1 

2(2i-l) , -2 ,4(i-l) 



i=0 

oo 



i=l 



2^ + l)^ 2 ^+*V 



2(1 - 1 2 ) 
1 + t 2 



1 + 



2t 2 



after some calculation. So the purity of the first excited mode is 



1-t 2 
1 + t 2 



1 + 



2t 2 



;i + t*y 



- n gs,l 



1 + 



2t 2 



2J1jJ\x 



;i + t 2 ) 2 J io + /2\ 2\u + fi< 

Since the other two modes contribute with two ground state factors, the total purity can be 
written as IIf s = rifs^ngs^IIgs^. For the purity parameter, one obtains finally 



Sera 1 



2\ 3 



and Sf s = 1 



1-f 
1+t 2 



2\ 3 



1 + 



2t 2 



'i + t*y< 



2\3 



= Sp 



2t 2 (l -t 2 ) 
(1+t 2 ) 5 ' 



1 + t 2 , 
In conclusion, Sf s < s gs . 

At long last, we may go back to Moshinsky's starting point, the assessment of electron 
correlation, only now for the excited state. The Hartree-Fock approximation for the relevant 
mode, in view of (8), is of the form 

WuF(r 1 ,r 2 ;pi,p 2 ) = ^ [W 00 (r 1 ;p 1 )W 11 (r 2 ;p2) - W 01 (r 1 ;p 1 )W w (r 2 \p 2 ) 

- W 10 (ri;pi)W 01 (r 2 ;p2) + W^ii(ri;pi)Woo(r 2 ;p2)], 
where W 00 (r;p) = -e'^-P 2 /^ W n (r;p) = -e'^-^^irjr 2 + p 2 /rj - \), 

7T 7T 

with their corresponding interferences. Remember that J Wij dr dp = 5^. In intracule-extra- 
cule coordinates: 

W HF {R,r;P,p) = ^{ r] r 2 +p 2 /r 1 - l) e -vR 2 -P 2 h-^-vVn. 

The parameter r] is determined by minimization. The mean value of the energy predicted 
by this function is: 

E ^ = \j (P 2 + oo 2 r 2 )[W 00 (r;p) + W n (r;p)} dr dp - 5 - j (n - r 2 f W KF (1, 2) dld2 

u 2 + 3/i 2 



35 

IT — V + 
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The minimum dE/drj = occurs when 77 = |\Al> 2 + 3/i 2 . Therefore, the energy predicted 
by Hartree-Fock is a/cj 2 + 3/i 2 . Thus, the "correlation energy" for the lowest excited state 
of harmonium is: 

E c , is = E & - E RF = ^i±^ _ ^ u 2 + 3/i 2 - 2 vV + /i 2 )/2 ~ - ^ • 

z 64 



Thus, the relative correlation energies are 




Figure 3: Relative correlation energy of the singlet and of the triplet excited mode. As 
expected, the relative correlation energy for the singlet is greater than for the triplet for 
small values of the coupling. At 5/u 2 ~ 0.67 the order is inverted. 

Figure 3 shows the exact dependence of the relative correlation energy for both systems 
as a function of 5/u 2 . The relative correlation energy for the singlet is greater than for the 
triplet, just as the purity parameter for the singlet is greater than the one for the triplet. 
At 5/u 2 = 0.67 the relation between these two quantities changes and the relative correlation 
energy for the triplet is greater than for the singlet. Note however that the entropy depends 
only the behavior of the occupation numbers, while the correlation energy has to do with 
the natural orbitals as well. Such a nice proportionality as (14) fails for the triplet state. 

Finally, Figure 4 shows the difference between the exact profile 1-density and the Hartree- 
Fock profile 1-density for the harmonium triplet, Phf(^) := / Wbf{Tj T%p\^p-i) dpi drv dp^. 
This description goes back to the Coulson-Neilson classic paper [46] on the helium Coulomb 
system. The "Moshinsky's hole" observed in the neighborhood of r = graphically shows 
the Hartree-Fock underestimation of the mean distance between the fermions, for the excited 
configuration of harmonium as well. 
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-0.2 



Figure 4: Moshinsky's hole for the triplet: (p(r) — pn F (r)) /co 1 ^ 2 as a function of u l l 2 r. 
9 Conclusion 

From the very beginning of quantum mechanics, the fundamental state of harmonium has 
provided a useful playground for learning about such questions as correlation energy, entan- 
glement or hole entropy (including black hole entropy). Here, for the first time, we rather 
exhaustively analyze the (spin triplet) first excited configuration of harmonium, particularly 
the behaviour of its occupation numbers and natural orbitals. This is a different chemical 
species altogether, due to the antisymmetric character of the orbital wave function. When 
exactly reconstructing d la Lowdin-Shull-Kutzelnigg the two-body density as a functional 
of the one-body density, instead of the sign dilemma (already solved by two of us) for the 
lowest-energy state, we find, as expected on general grounds, an ambiguity in the choice of 
natural orbitals. 

Also as expected, in the triplet case the first occupation number plays a more dominant 
role than for the singlet, up to fairly high values of the coupling parameter, t < 0.4. Thus, 
within this range, modeling the excited configuration as a Hartree-Fock state introduces a 
lower error than doing so for the ground state. In parallel, the linear entropy of the first 
excited configuration is lower than that of the ground state, and the relative correlation 
energy for the excited state stays below that of the ground state for such values of the 
coupling. The order reverses at higher values of t. 
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